Vibrational and Thermal Properties of ZnX (X=Se, Te): Density Functional Theory 

(LDA and GGA) versus Experiment 
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We calculated the phonon dispersion relations of ZnX (X=Se, Te) employing ab initio techniques. 
These relations have been used to evaluate the temperature dependence of the respective specific 
heats of crystals with varied isotopic compositions. These results have been compared with mea- 
surements performed on crystals down to 2 K. The calculated and measured data are generally in 
excellent agreement with each other. Trends in the phonon dispersion relations and the correspond- 
ing densities of states for the zinc chalcogenide series of zincblende-type materials are discussed. 

PACS numbers: 63.20.-e, 63.20.dk, 63.20.D-, 68.35.bg, 65.40.Ba, 71.55.Gs, 71.70.Ej 



I. INTRODUCTION 



In an ongoing effort to characterize the vibrational 
and thermal properties of multinary semiconductors and 
semimctals of current interest we have recently focused 
our attention on binary and ternary chalcogenides, such 
as HgX (X=S, Se, Te), CuGaS 2 or ZnS in the zinc-blende 
(zb) and rocksalt structured— Most of our past experi- 
mental work has been concerned with the specific heats 
in the low-temperature regime where deviations from the 
Debye-T 3 power law can be easily studied. Such devia- 
tions arise from low-energy phonon branches with low 
dispersion, typically zone-boundary transverse acoustic 
(TA) phonon branches. Ab initio calculations of the spe- 
cific heat (C„) based on the electronic structure derived 
from local density functionals have been carried out and 
compared in detail with the experimental results. They 
were found to account rather well for the experimental 
data and, in some cases, confirm our experimental re- 
sults where large deviations from previously published 
measurements were encountered, e. g. for GaN£ For 
most of our investigated systems, by selective isotope 
substitution we could effect incremental variations of the 
low-temperature heat capacities^ The small changes of 
the heat capacities induced by selective isotope substi- 
tution arc generally very well described by density func- 
tional theory (DFT) and were advantageously used to 
contrast the calculations with experimental data. The 
ab initio calculations were done employing the two most 
common approximations for the density functional of 
the exchange and correlation energy, the local density 
approximation (LDA) and the generalized gradient ap- 
proximation (GGA) implemented in various DFT codes. 
Usually, LDA and GGA results bracket the experimen- 
tal data. For example, LDA tends to underestimate the 
lattice parameters and the cell volume, while the GGA 
approximation overestimates them. This has lately been 
discussed by us, for example, for ZnS, the cubic mer- 
cury chalcogenides HgX (X=S, Se, Te) and cinnabar (a- 
HgS). For systems including heavy elements we have also 
tested the effect of including spin-orbit (SO) interaction 



into the DFT calculations. For a-HgS, SO coupling was 
found to give a negligible effect on the thermal prop- 
erties. However, the heat capacity and also the phonon 
dispersion relations of the semimetal Bi can be noticeably 
better described if SO interaction is included in the DFT 
calculations £ SO interaction effects on the thermal prop- 
erties are less pronounced for Sb, another semimctalfS 
The relativistic electronic structures which results from 
our calculations enabled us to study a number of interest- 
ing band structure effects such as linear SO splittings of 
the valence band near the origin of the Brillouin zonei* 3 - 
In addition to this program, also various DFT implemen- 
tations (ABINIT and VASP) have been compared with 
each other and probed against the experimental data. 

In the present contribution we extend our investiga- 
tions to the thermal and vibrational properties of the 
binary semiconductor ZnSe and ZnTc. ZnSc and ZnTe 
crystallize at normal pressure and temperature in the cu- 
bic zb structure. ZnSe is a wide-band-gap II- VI semicon- 
ductor with a gap of 2.8 eV while the band-gap of ZnTe 
is smaller and amounts to 2.2 eV. Because of its supe- 
rior wide range infrared low-temperature optical trans- 
mittancc ZnSc has gained technological importance for 
application as infrared transparent optical components, 
such as windows. ZnTe can be easily doped and found ap- 
plications in optoelectronics, e.g. for blue light-emitting 
diodes and solar cells. 

The vibrational and thermal properties of ZnSe and 
ZnTe have been the subject of a number of studies. Ini- 
tially the comparison of experimental data such as the 
phonon dispersion branches of ZnSe and ZnTe with in- 
elastic neutron scattering (INS) (Ref. [ToUXlT ) or the heat 
capacities with theoretical approaches remained on the 
semiempirical leveled— These early investigations have 
been extended recently to ab initio LDA and GGA cal- 
culations using up-to-date implementations of the DFT 
method. LDA results have been compared with the INS 
phonon energies and the heat capacities and found to be 
in good agreement.— ~— Hamdi et al. have also calculated 
the structural parameters of ZnSe and compared them 
with experimental and other available theoretical dataJ^ 
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Recently, Tan et al. have used DFT-GGA calculations 
to study the pressure induced phase transitions from the 
zb to the rocksalt structure.— They also performed cal- 
culations of the phonon dispersion relations and heat ca- 
pacities of ZnSe and ZnTe with the DFT-GGA CASTEP 
code.— Unfortunately, the results of these calculations 
have not been compared with the available experimental 
data in more detail. 

Investigations of isotope substitution of the atomic 
masses of either Zn/(Se,Te) or simultaneously of both 
constituents simultaneously, on the properties of ZnX (X 
= Se, Te) appear to be rather scarce. To the best of 
our knowledge, the only available studies are a theoret- 
ical investigation of the temperature and mass depen- 
dence of the lattice constants, the temperature depen- 
dence of the linear thermal expansion coefficient and the 
mode Griineisen parameters upon isotopic substitution 
in ZnSe using perturbation theory in a density-functional 
framework.— Also, available is an experimental investi- 
gation of the Raman and excitonic spectra with different 
isotopic compositions.— In Ref. Hi] an interesting depen- 
dence of the linewidth of the Raman spectra of ZnSe on 
isotope disorder was discovered. 

We have recently demonstrated that changes of the 
specific heats due to the isotope substitution, even 
though sometimes very small, can be revealed by 
careful experiments even on very small crystalline 
samples, 3 Especially, we have shown that the logarith- 
mic derivatives with respect to the atomic masses m i; 
din (C p /T 3 )/dln rrn are a very sensitive and useful tool 
to highlight small alterations of the specific heats since 
the logarithmic derivative eliminates, to a large extent, 
systematic experimental errors. In binary and multinary 
compounds the logarithmic derivatives particularly al- 
lowed us to reveal the influence of variations of the atomic 
masses on the phonon spectrum and to compare with the 
calculated phonon density of states (PDOS). We have 
also proposed and tested a sum rule which relates the 
sum of the logarithmic derivatives of C p /T 3 with respect 
to the atomic masses with the logarithmic derivative with 
respect to temperature^ 

The low-temperature heat capacities of ZnSe have been 
the subject of two preceding experimental studies car- 
ried out in the 1970s by Irwin and La Combed and by 
Birch<2i However, the two previous heat capacity data 
sets just overlap in the temperature range of the C p /T 3 
maximum which for ZnSe occurs at about 20 K: in order 
to enable a meaningful comparison, e.g. with the loga- 
rithmic derivatives, and to test the sum rule a repetition 
of the measurements on a ZnSe sample with natural iso- 
tope composition was found necessary. We therefore car- 
ried out measurements of the heat capacities on a series 
of samples with various isotope combinations with special 
emphasis on the low-temperature regime. For ZnTe heat 
capacity data are available only down to ~15 K . 23 ' 25 ' 26 
The data from the different authors deviate significantly 
from each other, especially in the temperature regime 
where the maximum in Cp/T 3 is expected. We therefore 



repeated and extended the heat capacity measurement on 
ZnTe down to ~ 2 K. Our measurements reveal clearly 
the maximum in C p /T 3 at ~14 K. Additionally, our low 
temperature heat capacity data enable us to extract very 
reliably the Debye temperatures for T — > K, ODebye(O) 
for ZnSe and ZnTe. 

The comparison of our experimental data and litera- 
ture data of the phonon dispersion obtained by INS ex- 
periments with the calculations reveal significantly better 
agreement with the LDA while the GGA data deviate no- 
ticeably from the experiments. The analysis of the total 
and the partial phonon densities of states and the com- 
parison with that of ZnS (Ref. 0) highlights the close 
similarities of the PDOS of ZnX (X = S, Se, Te). It also 
reveals the compression of the frequency scales with in- 
creasing masses of the chalcogen atom. The PDOSs pro- 
jected on the Zn and the X atoms shows an interesting 
trend on going from S to Te. While in ZnS the acoustical 
and optical phonons can be clearly attributed to vibra- 
tions of the Zn and S atoms, respectively, the situation 
is reversed for ZnTe. Here the optical branches emerge 
essentially from vibrations of the lighter Zn atoms. For 
ZnSe, with almost equal masses of the Zn and the Se 
atoms, the acoustical and optical phonon branches arise 
from combined vibrations of Zn and Se, respectively. 



II. THEORETICAL DETAILS 

For the calculations reported here we have made use 
of density functional theor y 27 ' 28 as implemented in the 
ABINIT packago 2 ^. Two different exchange-correlation 
functionals (LDA, GGA) were tested but most of our 
results arc based on the LDA approach, basically because 
it gives a much better agreement with experimental data. 

In the ABINIT package we have utilized a linear re- 
sponse approach 3 -^— together with an iterative min- 
imization norm-conserving pseudopotcntial plane-wave 
method^ These pscudopotentials are single projec- 
tor, ordinary norm conserving, based on the Troullier- 
Martins method^ - — 12 and 6 "valence" electrons were 
used for Zn and Se, respectively. A 40 Ry cutoff was set 
for the plane wave expansion and an 8 x 8 x 8 k-point 
grid regularly shifted along four different directions, with 
a total of 60 k-points for the ground state calculation 
and 2048 k-points for the vibrational response part. The 
calculations converged well: increasing the plane-wave 
cutoff to 50 Ry and the k-points mesh to 10 x 10 x 10 
yielded, on average, a change of less then two percent 
in the phonon frequencies. For the phonon frequencies 
we have used a q-mesh of 12 x 12 x 12, to guarantee a 
good coverage of the dispersion relations. The dynami- 
cal matrices were obtained from perturbation theor y 31 ' 32 
and a Fourier interpolation was employed in order to in- 
crease the mesh sampling and to improve the description 
of quantities such as vibrational density of states and 
heat capacities. 
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III. EXPERIMENTAL 

zb-ZnSe and zb-ZnTe crystals with different isotopic 
compositions: nat Zn nat Sc/Tc i.e., Zn and Se/Te with 
natural isotope composition, 64 Zn 76 Se, 64 Zn 80 Se, 68 Zn 
76 Se, and 68 Zn 80 Se, 68 Zn nat Te, and nat Zn 130 Te, were 
grown by vapor phase transport as described in detail 
previously.—. Larger zb-ZnTe crystals with the natural 
isotope composition were obtained by a modified Pipcr- 
Polich method in a semi-open quartz-glass ampoule i 38 ! 39 
The heat capacities of crystalline pieces of these samples, 
typically of ~20 - 100 mg weight, were measured between 
2 and 280 K with a physical property measurement sys- 
tem (Quantum Design, San Diego, CA) as described in 
detail in Ref42 



IV. LATTICE DYNAMICS 

In Figure [Ha) and (b) and we display the calculated 
phonon dispersion relations of ZnSe along selected di- 
rections of the Brillouin zone (BZ). For comparison we 
have plotted our LDA and GGA results and the INS data 
reported by Hennion et al— The LDA calculations are 
overall in good agreement with the experimental findings, 
as has also been observed by Dal Corso et al— The re- 
sults of the GGA calculation deviate markedly, especially 
for the optical branches. 

Figure HJa) and (b) show the comparison of the LDA 
and GGA calculations for ZnTe with the INS data col- 
lected by Vagclatos et al.— As already observed by Dal 
Corso et al. the agreement of the experiment with the 
DFT-LDA calculations is excellent, especially for the 
acoustical branches. The longitudinal optical branches 
are equally well reproduced by the calculations, while the 
transverse optical branches fall short by about 5%. The 
DFT-GGA results agree well with the optical branches 
while the energies of transverse acoustical branches lie 
systematically below the experimental data. 

The total PDOS and the PDOS projected on the two 
constituting elements as obtained by the LDA calcula- 
tions, are displayed in Fig. [3] For comparison we also 
show the PDOS of ZnSe and ZnTe derived by Talwar et 
al. from a rigid- ion-model^ In Figure 0] we also added 
the PDOS of ZnS as published by us recently! 4 - Very 
similar to the PDOS of ZnS the PDOSs of ZnSe and 
ZnTe exhibit two dominating features, one between 50 
and 100 cm -1 , arising essentially from acoustic phonon 
branches, and a sharp double peak at higher energies cor- 
responding to the the optical phonons. The energies of 
the latter peaks decrease markedly as one substitutes S 
by Se and Te. These two main peak groups are separated 
by a gap in which a small mid-gap peak is located. While 
the low-energy peak in the PDOS of ZnSe is largely in- 
dependent of the computational approach the mid-gap 
feature and the sharp optical phonon double-peak ob- 
tained from the GGA calculations are down-shifted, the 
latter by about 10% as compared to the LDA result. Ex- 
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FIG. 1: Phonon dispersion relations of ZnSe as obtained by 
DFT-LDA (a) and DFT-GGA (b) calculations. The circles 
represent the inelastic neutron scattering data measured by 
Hennion et al.— 



ccpt for the mid-gap peak the LDA PDOS is generally in 
good agreement with results of the rigid-ion-model used 
by Talwar et al.— 

It is of interest to compare the projected PDOS of the 
three systems ZnX (X =S, Se, Te) (see Fig. i}. For the 
ZnS the acoustical and the high energy optical phonon 
branches arise essentially from vibrations of Zn and S, 
respectively. Due to the reversal of the masses of the two 
constituents in ZnTe (Zn: ~65 a.m.u.; Te: ^128 a.m.u) 
the situation is reversed for ZnTe. There the acoustical 
phonons arise from Te vibrations and the optical ones 
from vibrations of the lighter Zn atoms. Due to the very 
similar atomic masses of Zn and Se acoustical and optical 
phonon branches for ZnSe comprise mixed Zn and Se 
vibrations. 



V. HEAT CAPACITY 

By using the calculated PDOS(w) we obtained the free 
energy, F(T) and the heat capacity (at constant volume), 
C V (T) according to— 
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FIG. 2: Phonon dispersion relations of ZnTe as obtained by 
DFT-LDA (a) and DFT-GGA (b) calculations. The circles 
represent the inelastic neutron scattering data measured by 
Vagelatos et aZ..— 



FIG. 3: One-phonon density of states (PDOS) of (a) ZnSe and 
(b) ZnTe as obtained by DFT-LDA and DFT-GGA calcula- 
tions. The diamonds represent data by Talwar et aL.— The 
inset displays the high energy section in an enlarged scale. 



F(T) = - J |— + k B T]n[2n B {u)}jp{u))duj, (1) 
o 

a(T) = -r(g) v , (2) 

where fcs is the Boltzmann constant and tib the Bosc- 
Einstein factor. The cut-off in PDOS(cj) at the highest 
phonon frequency, defines the upper limit of integration. 

Figure [5] and [6] display the calculated heat capacities, 
C p ~ C v (T) of ZnSe and ZnTe for different isotope mass 
compositions in the C p /T 3 (T) representation which is 
chosen to emphasize the broad peak at low tempera- 
ture. For ZnSe the maximum is found at 17.2 K and for 
ZnTe it moves down to 13.5 K. The broad peak in the 
C P /T 3 (T) representation of the heat capacity has been 
demonstrated to arise from regions of low dispersion of 
zone-boundary acoustical phonons which also give rise to 
the low-energy peaks in the PDOS.— They appear at 68.6 
cm -1 and 55.1 cm -1 for ZnSe and ZnTe, respectively. We 
again find a ratio ~ 6 for the energies of the acoustical 
phonon peaks and the peak temperature in the C P /T 3 {T) 
representation, close to what has been typically observed 
for several semiconductors before.— 

The peak value of C p /T 3 (~ 13.5A") for ZnTe is by al- 



most a factor of two larger than that of ZnSe at 17.2 
K, reflecting the lattice softening when Se is replaced 
by Te. Isotope mass increases lead to small growth 
of the heat capacities, which become clearly visible as 
slight changes of the maximum value of the C P /T 3 (T) 
peak. Below we will discuss the logarithmic derivatives of 
C P /T 3 (T) which reveal that the variation of the isotope 
masses changes the heat capacities in different tempera- 
ture regimes. Figure [7] displays the logarithmic deriva- 
tives, din (Cp/T 3 ) j 'din nii, with respect to the isotope 
mass, of the two constituting elements, Zn and either 
Se or Te, respectively. For completeness we have also in- 
cluded the analog graphs for ZnS<^ For the latter there 
is a clear differentiation between the derivative vs the 
masses of Zn and S. Similar sharp peaks at low temper- 
ature, with that of S significantly smaller in magnitude 
than that of the Zn derivative, have been attributed to 
a considerable S component in the TA vibrations which 
have mainly Zn character. The characteristic S contri- 
bution due to high energy S vibrations shows up as a 
broad peak centered near 100 Similar derivatives are 
obtained for ZnSe and ZnTe. However, depending on the 
relative masses of the constituents Zn and Se/Te the mag- 
nitude and the maximum temperatures of the peaks shift. 
For ZnSe the low-temperature peaks in the derivative are 
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FIG. 4: Comparison of the one-phonon density of states 
(PDOS) of (a) ZnS (Ref. |), (b) ZnSe and (c) ZnTe as ob- 
tained by DFT-LDA calculations. The decomposition into 
contributions from Zn and from the chalcogen is also given.— 



of equal magnitude and a pronounced high-temperature 
peak is missing. This finding reflects the mixing of Zn- 
and Sc-related vibrations due to the rather similar atomic 
masses of Zn and Se. For ZnTe the situation is reversed. 
Here the pronounced peak at low temperature is associ- 
ated to Te vibrations, while the vibrations of the lighter 
Zn atoms give rise to a shallow broad peak centered at 
~ 70 K. Again, as in ZnS, there is some hybridization of 
Zn with Te vibrations and a noticeable Zn component in 
the TA vibrations which now have mainly Te character. 

In previous work we have demonstrated that for low 
temperatures, T — > 0, the logarithmic derivatives are re- 
lated to the ratios of the atomic mass to the molar mass 
according to^£^ 
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FIG. 5: Calculated low-temperature heat capacity of ZnSe 
(ABINIT DFT-LDA) assuming different masses for the Zn 
and the Se isotopes, as given in the inserts (increasing isotope 
mass from bottom to top). 
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(3) 



m x 



2 m Zn + m x 



0.82(Se);= 0.99(Te). 



(4) 

The dashed lines at low temperatures in Figure [7] rep- 
resent a free-hand extrapolation for T — >• so as to agree 
with these values. 

Our low temperature heat capacities of ZnSe and ZnTe 
collected on samples with the natural isotope composi- 
tion are summarized and compared with the results of the 
DFT-LDA and GGA calculations in FigureH The C p /T 3 
representation reveals low-temperature peaks at 17.0 K 
and at 13.5 K, for ZnSe and ZnTe, respectively. For ZnSe 
the data by Irwin et al£2- and Birch^l are in good agree- 
ment with our more complete data. Especially, for ZnTe 
we could resolve the diverging low-temperature heat ca- 
pacity results obtained before by Gavrichev et al££- and 
Irwin et al^ and extend the temperature range down to 
2 K. This was found necessary in order to clearly establish 
the position and magnitude of the C p /T 3 maximum and 
to enable the calculation of the logarithmic derivatives, 
discussed below. 
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FIG. 6: Calculated low-temperature heat capacity of ZnTe 
(ABINIT DFT-LDA) assuming different masses for the Zn 
and the Te isotopes, as given in the inserts (increasing isotope 
mass from bottom to top). 



The Debye temperatures, 0£> e t ye (0), for T — s- K are 
obtained from the experimental data by fitting a power 
law 



C p /T = (3T 2 + ST 4 

to the low-temperature range 2 K < T < 10 K. 
For ZnSc we obtained 



(5) 



0Debye(O) = 289 (2)K for ZnSe, 

which is by about 6% larger than the value reported 
by Birch^i and somewhat closer to the result (278.5 K) 
derived from the elastic constants by Lee42 

For ZnTe we found 



0Dcbyc(O) = 230(2)K for ZnTe 

which is close to the value from elastic constants measure- 
ments 225.3 but significantly larger than the values 
estimated from x-ray measurements and quoted in review 
articles (182 K - 210 K)M^1 

We finally note that the ratio of the Debye- 
temperatures, 1.26, is very close to the square-root of 
the inverse atomic mass ratio of the chalcogenide atoms, 
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FIG. 7: (Color online) Logarithmic derivatives vs the masses 
of the constituents, Zn (solid red) and X (X = S, Se, Te; 
from top to bottom, see inset) (solid blue) calculated from 
the LDA heat capacity results displayed in Figs. [5] and [6] 
The ZnS data (a) were taken from Ref. Q. The dashed part 
of the calculated curves have been extrapolated by free hand 
so as to agree with the value of Eqs. © and @ at T = 0. 



As already observed for the phonon dispersion, the 
agrement of the heat capacities of ZnSe and ZnTe with 
the LDA calculations is better than with the GGA re- 
sults. Especially in the temperature regime of the C p /T 3 
maximum overestimates/underestimates the experimen- 
tal data by about 20%/50% for ZnSe and ZnTe, respec- 
tively. 

The logarithmic derivatives vs the mass of the isotopes 
of Zn and Se are summarized in Figure [5] Those for Zn 
and Te are displayed in Figure I10I The agreement with 
the LDA calculations in magnitude and position of the 
low temperature peaks is very satisfactory. For ZnTe, the 
logarithmic derivative with respect to the isotope mass of 
the Zn atoms clearly exhibits the broad maximum cen- 
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FIG. 8: Our experimental results and literature data of the 
low-temperature heat capacities of (a) ZnSe and (b) ZnTe 
compared with the ABINIT DFT-LDA and GGA results, 
(black) solid and (red) dashed lines, respectively. For the 
references of the literature data see text. 



FIG. 9: Logarithmic derivative of C p /T 3 of ZnSe vs the iso- 
tope masses of (a) Zn and (b) Se. (o) experimental results, 
(red) solid line LDA calculations. The (red) dashed lines are 
free-hand extrapolations (see Fig. . 



tered at ^50 K due to phonons involving low-energy Zn 
vibrations. 

In previous work, on monatomic crystals we derived a 
connection between the logarithmic derivative of C p /T 3 
versus T and the corresponding derivative vs the isotopic 
massif For binary and ternary materials, a similar con- 
nection was shown to hold provided one adds the two or 
three logarithmic derivatives with respect to each of the 
isotope masses.^22 

Accordingly, for the relation of the temperature and 
isotope mass derivatives in a two component system we 
obtain the expression^ 



1 dln(C p (T)/T 3 ) = * dHC p (T)/T 3 ) 

2 [ rflnT ' ^ dlrnrii 

i—l 



where rrii (i =1, 2) are the masses of the two constitut- 
ing elements, i.e. Zn and X = Se and Te, respectively. As 
demonstrated in Figure [TTJ the relation given in Eq. © 
is nicely fulfilled for ZnSe and ZnTe, thus lending further 
support for the proposed sum rule. 



VI. CONCLUSIONS 

Up-to-date ab initio electronic structure calculations 
have now matured and became a powerful tool to study 
the thermodynamic properties of crystals. Here we 
demonstrate this for the binary zinc chalcogenidcs, ZnSe 
and ZnTe. Both materials are technologically important 
semiconductors used for optical applications. We have 
used DFT-LDA and GGA codes to calculate the phonon 
dispersions and the heat capacities down to lowest tem- 
peratures. In the experimental part we extend the avail- 
able experimental basis and comprehensively determine 
the heat capacities of ZnSe and ZnTe over the tempera- 
ture range from 2 K to room temperature and have com- 
pared the results with the calculations. Generally, better 
agreement of the experimental data is found with the 
DFT-LDA results. Additionally, we measured the heat 
capacities of isotope substituted samples in order to vary 
the atomic masses of Zn and Se and Te and to probe 
the effect of selective isotope substitution on the phonon 
spectra and the thermodynamic properties. These data 
are compared with the results of the calculations and 
found to be in very good agreement. A comparison of 
ZnS, ZnSe and ZnTe instructively reveals the variation 
of the phonon spectrum with the variation of the atom 
mass of chalcogcnidc atoms. 
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FIG. 10: Logarithmic derivative of C p /T 3 of ZnTe vs the 
isotope masses of (a) Zn and (b) Te. (o) experimental results, 
(red) solid line LDA calculations. The (red) dashed lines are 
free-hand extrapolations (see Fig. . 
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